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Abstract. - The Lassig- Wiese (LW) field theory for the freezing transition of random RNA 
secondary structures is generalized to the situation of an external force. We find a second-order 
phase transition at a critical applied force / = f c . For / < f c forces are irrelevant. For / > f c , 
the extension £ as a function of pulling force / scales as £(/) ~ (/— /c) 1 ^ 7-1 - The exponent 7 is 
calculated in an e-expansion: At 1-loop order 7 = e/2 = 1/2, equivalent to the disorder-free case. 
2-loop results yielding 7 = 0.6 are briefly mentioned. Using a locking argument, we speculate that 
this result extends to the strong-disorder phase. 



Introduction. — RNA is a heteropolymer con- 
^^structed from four different nucleotides A, C, G and U 
""^located on a sugar-phosphate polymer backbone, fn so- 
. — lution, a single RNA strand bends back onto itself and 
refolds into a configuration of loops, stems and termi- 
QTtiating bonds, due to formation of Watson-Crick pairs 
^A-U and C-G from bases located on different parts of 
• — the strand. Together with environmental conditions like 
temperature and ionic concentration, the primary struc- 
H ip,re (base sequence) determines the most probable base- 



pairings, known as secondary structure, which then de- 
termines the most probable spatial conformation (tertiary 
structure) [1,2]. Unlike protein folding, which exhibits 
a strong interdependence between secondary and tertiary 
structure [3], RNA folding may be studied at the level of 
secondary structures due to a clear separation of energy 
scales. 

Since the pioneering work of Bundschuh and Hwa [4], 
several authors have studied the statistical physics of RNA 
secondary structures for random sequences [5-8]. It is 
commonly believed that these systems undergo a freez- 
ing transition upon lowering the temperature. Based on a 
replica approach, Lassig and Wiese [9], and David and 
Wiese [10] have recently developed a systematic field- 
theory formulation for this phase transition in terms of 
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interacting random walks (RW). The critical exponents 
characterizing pairing statistics and replica overlap were 
computed within a 2-loop renormalization analysis, and 
found to be remarkably close to numerical simulations 
[4,11]. 

An interesting way to probe RNA chains is to study its 
behavior under an external pulling force (see fig. la for an 
illustration). Recently, force-extension curves have been 
measured, by attaching beads to the RNA-molecule and 
pulling on it, using an optical trap [12, 13]. For homopoly- 
mers, the competition between structure formation and 
denaturation of the RNA strand leads to a second-order 
phase transition at a critical force / = f c [11]. For / < f c 
the strand is still in a collapsed phase, while for / > f c 
it is in a extended "necklace phase" with a macroscopic 
extension (end-to-end distance). While quite some litera- 
ture exists about the subject (see e.g. [14, 15]) there is of 
today no theory to compute the characteristic exponent 7 
of the force-extension characteristics for disordered RNA 
strands at the transition. In this letter we fill this gap. We 
propose an extension of the RW field theory for random 
RNA by including an external pulling force /. We show 
that within our theory the force is renormalized by the 
quenched disorder, hence the exponent 7 is modified with 
respect to its mean-field value 70 = 1/2. Conversely, we 
argue that the disorder coupling is not renormalized by 
the force at any order in perturbation theory. We use per- 
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Fig. 1: (a) Illustration of the planar structure of RNA under 
tension, (b) rainbow diagram and (c) corresponding height re- 
lief. 



turbative renormalization of the new theory to compute 7 
to 1-loop order. Finally, we comment on 2-loop results. 

The model. — We consider an RNA strand with L 
bases labeled with indices s = 1, . . . , L. Similarly, we use 
the index s to label a backbone segment between adja- 
cent bases s and s + 1. A secondary structure S is a 
set of base pairs (s,t), s < t. We only retain so-called 
planar structures S: any two different base pairs (s,t), 
(s',t') are either independent s < t < s' < t' or nested 
s < s' < t' < t. S can be represented by a diagram of 
arches (fig. lb). Planarity implies that the arches do not 
cross each other. We introduce the contact operator $ 
defined by $(s,t) = 1 if (s,i) 6 S, and otherwise [9]. 
Moreover we define the height field h on the segment r by 
M r ) := E s <r St> r ^( s ' coun ting the number of arches 
over r. This leads to an identification of each open pla- 
nar secondary structure S with a height function subject 
to boundary conditions h(0) = h(L) — 0. A segment be- 
tween the bases r and r + 1 belongs to the free part of the 
structure if and only if h(r) = (fig. lc). 

In order to develop a statistical mechanics model, we 
have to assign to each structure S an energy E[S]. We 
assume that it may be written as a sum of the contribu- 
tions from the formation of base pairs -E pa ir[<S'] and from 
the external force Ef orce [S]. Bond formation between the 
bases at s and t involves a pairing energy rj(s,t) which 
in general depends on the nature of the pairing partners. 
We sum over all pairing energies of base pairs in S and ob- 
tain E pail [S] = £ s<t r?(M)$(s,i) [4]. The energy due to 
the external force, Ef orcc [S], depends on the spatial con- 
figuration of the free part of the strand and its elasticity. 



We assume that every free backbone segment aligns with 
the force, hence the energy is proportional to the force 
/ times the number of monomers in the free strand [6]. 
Thus, we neglect any elasticity and entropic effects for 
the free segments and for the bonds which terminate loop 
structures (figure la). By analogy with the contact opera- 
tor <!>(s,t), we introduce a free-strand operator A(r) such 
that A(r) = 1 if h(r) — 0, and otherwise. This allows to 
write £ forco [S] = -/£ r A(r) . 

Having defined the energy of a given secondary struc- 
ture, we proceed to study the partition function 



z = £ ex p 

ses(L) 



•5>(M)*(M) + /£ A ( S ) 



(1) 



where S(L) denotes the set of all possible planar sec- 
ondary structures with L bases. Before considering ran- 
dom RNA chains, we briefly review the properties of the 
partition function in the case of uniform pairing ener- 
gies rj(s,t) = ijo (that we may take = 0). For / = 
one deduces from the height picture that the problem is 
equivalent to the statistics of a RW on the positive real 
axis h > 0. The partition function is Zl oc L~ Pa with 
Po = 3/2 the characteristic exponent of first return. This 
leads to a pairing probability for the base pair (s, t) scaling 
like p(s,i) = ($(s,t)) « [\t - s\(L - \t - s|)]~ Po . Switch- 
ing on the force / > amounts to adding an attractive 
short-range potential at the origin h = 0. This is a well- 
known problem of statistical mechanics. For instance it 
describes surface wetting transitions in 1 + 1 dimensions 
(see e.g. [16]). For forces / larger than a critical force f c , 
the RW is bound to the origin h = whereas for f < f c 
it is unbound (i.e. free to wander far away from h = 0). 

In fact, the / = problem can be mapped onto a free 
RW r(s) in d = 2po = 3 dimensions; the height field is the 
modulus h(s) = \r(s)\. The inclusion of the short-ranged 
attraction in the case / > corresponds to a short-ranged 
attractive potential at r = 0. In the continuum limit, the 
action of this model reads 



S 3D {r(s)}=£d S Q[r( S )] 2 



-fS 3 (r(s))) , (2) 



and describes the pinning of a RW by an attractive impu- 
rity at the origin. 

We now exploit this analogy to extend our analysis to 
random RNA structures. Numerical simulations [4] sug- 
gest to model sequence disorder by independent Gaussian 
random binding energies rj(s,t), 



rj(s,t)=T)o, r)(s,t)r)(u,v) -r] = a S(s-u)S(t-v) (3) 

Following [9,10], we construct a field theory in the con- 
tinuum limit L — * 00. We perform a perturbative ex- 
pansion in the disorder amplitude g — a 2 > 0, and the 
force strength /. To model disorder, we use the replica 
trick. Each replica a — l,...,n is represented by a 
RW r a (s) in an embedding space M. d (with dimension 
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d = 2p = 3). In fact, the explicit form of the pairing 
probability for uniform RNA suggests to consider closed 
RWs r a (s). Nevertheless, we can (and shall) use open 
RWs because they have been proven to lie in the same 
universality class and considerably simplify the calcula- 
tions [10]. Within the RW representation, the contact 
operator reads $ a (s,t) — S d (r a (s) — r a (t)). The aver- 
age over the disorder r\ generates an attractive interac- 
tion between the replicas. It is described by the overlap 
operator ty al 3(s,t) = $ a (s,i)$p(s,t) (counting the com- 
mon arches of the replicas a and j3 in the original pic- 
ture). By analogy with (2) we represent the operator A 
as A Q (s) := S d (r a (s)). The resulting action in the RW 
picture is 

S[{r a }]=Y, f\[*°(s)] 2 

a J 

-9Y.il -/E / A «^> w 

a <(3 JJs<t a=l J 

and generalizes the model of [10] (where / = 0). Before 
using perturbation theory, we generalize the model to di- 
mensions 2 < d < 3 [9, 10]. Setting e = d — 2 we find the 
canonical scaling dimensions dim g = e and dim / = e/2. 
The original theory corresponds to e = 1 . The generalized 
model is renormalizable at e = as the / = model [17]. 

Perturbation theory. — We represent the perturba- 
tive expansion of Z in /, g in terms of Feynman diagrams. 
The ^-vertex (disorder interaction vertex) is denoted 
by a double arch between a pair of replicas [9, 10] and the 
/-vertex A a (force interaction at r = 0) is depicted by a 
force insertion on a single replica: 



A Q 



(5) 



Non-planar diagrams involving the disorder interaction, 
see fig. 2(a), may be eliminated by introducing nxN pairs 
of auxiliary fields 7 and 7 in the action and by taking the 
large- N limit [10]. Since the force term only acts on the 
free part of the RNA strand, also diagrams of the type 
given on fig. 2(b) must be excluded. This is achieved by 
a similar "planarity constraint" that can be implemented 
using the same auxiliary fields. 



4> !-> — 
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Fig. 2: Non-planar diagrams. 

We now consider the partition function for a single (bun- 
dle of replica) RW with fixed end-points, or rather its 



Fourier transform Z({q' a } 7 {q^}; /, g), defined as 
Z({q' a ,q'>};f,g) 



[v[{r a }}e- s ^Y[V q l(0)V q l(L) 



(6) 



where the vertex operator V q (s) = exp(iqr a (s)) in- 
jects incoming external momenta {q' a } 7 {<?"} to each end- 
point of the replica a of the RW. The (regularized) di- 
mcnsionlcss integration measure is given by X^r^] = 
Hi d d r Q ^/(47ra 2 /m) d / 2 where a is an ultraviolet cut-off 
and m — 1/2 the mass of the Brownian particle asso- 
ciated with the RW. We further simplify the model by 
setting q' a = q', q'a = q" without loss of generality. 

The perturbative expansion in g and / leads to a 
systematic diagrammatic representation of Z. The dia- 
grams can be classified according to the number of repli- 
cas with force insertions. The absence of force inser- 
tions on the replica a leads to a momentum conservation 
3 d (q'a + Qa) (translation mvariance r a (s) -> r a (s) + r). 
Thus, the set of all possible diagrams is classified accord- 
ing to the number k of its external momentum conserva- 
tions: k — n conservations correspond to a diagram of 
the force-free theory, k = n — 1 conservations to a single 
replica subject to force insertions, etc. We group all dia- 
grams with k conservations into a restricted partition func- 
tion Z n _ k (q',q";f,g), so that formally Z(q',q";f,g) = 
X^<K<?' + q") k Z n -k(q',q";f,g)- S describes the force- 
free theory, Si the sector where the attractive short-range 
potential only acts upon a single replica, E 2 the sector 
where it acts upon two replicas, etc. In the following we 
focus on Ei(g, — q), since this simplifies the calculations. 
Perturbative expansion up to order two in / or g yields 

5i(q, -q;f,g) = 

fn : + f 2 n ! ! 



fg 



n(n — 1) ( / 1 /"V 7 + . 



+ 



n(n - l)(n - 2) 



+ •■ 



(7) 



There are four topologically different Feynman diagrams. 
The first contribution is 



/ 



(Ai(u) J] V?(0)V? q (L))o 

0<m<L 

Le- q2Ln , 



a=l 



(8) 



where ( }o denotes the average w.r.t. the free RW action 
So[r] (/ = .g = 0), with a proper subtraction of the (n — 
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1) translational zero- modes. The second diagram is UV 
divergent at e = 0. Isolation of the corresponding pole 
yields 



J J (Ai(u)Ai(u) JJ Vq(0)V" q (L)) 



0<u<v<L 



a=l 



\v — u\ ' e 



a<u<v<L 

= -— e - q2Ln + o(i) 



(9) 



The divergence comes from the short-distance behavior 
of the product of two A operators A(u)A(v) = \u — 

u — >v 

v\~ d / 2 A(v) + • • • . The third contribution is also UV di- 
vergent, 



(Ai(u)*i2(u,«;) J] Vf(0)V? q (L))o 



0<u<v<w<L 



a=l 



e -2q 2 (L-w+v), w _ v \-d £ -q 2 L (n-2) 



0<u<v<w<L 



= q --(Lq 2 - 1) e- q ' Ln + 0(1) 
as well as the fourth 



(10) 



0<u<L 
0<v<w<L 



a=l 



e -2q 2 {L-w+v)\ w _ v \-d e -q z L(n-2) 



0<u<L 
Q<v<w<L 

= q ~ (2Lq 2 - 1) e- q2Ln + 0(1) . 



(11) 



Renormalization. — We remove the UV divergences 
in the expansion (7) by formally taking po a s analytical 
regularization parameter. In order to eliminate the simple 
poles in e = 2p — 2 at e = 0, we define the renormalized 
theory through the renormalized action 

S " = E / ~fi - 9R^ a E // *^( s > *) 

- f R »- e/2 Z f [ A a (s) + 2nZi . (12) 

J s 



Here Z, Z g and Z / denote the wave- function, the coupling 
constant and the force counterterms respectively. Zi ac- 
counts for boundary effects since we deal with open RWs 
[10]. The coefficients of their development in fx, gn con- 
tain the leading poles in 1/e. Furthermore, we have intro- 
duced the renormalization mass scale \i ~ 1/L. From di- 
mensional analysis we deduce the relations between renor- 
malized and bare r, q, g and /: For the field r = Z 1 / 2 r^, 
q = Z~ 1//2 g#, for the couplings g — gF>p,~ e Z g Z~ d and 
/ = /flA i_e/ ' 2 ^/^ _d ^ 2 - The renormalized partition func- 
tions are related to the bare ones via 



{Qr,Qr,JR,9r) 



= z- kd ' 2 ( 



-2nZ 



^E k (q',q",f,g) . (13) 



The prefactor on the r.h.s. takes into account the (n — k) 
zero modes, boundary effects and the change of normaliza- 
tion in the integration measure. Consistency of the theory 
requires cancelation of all divergences upon renormaliza- 
tion for each individual Efc. In particular, renormaliza- 
tion of the force-free terms So has been performed previ- 
ously [10] and yields the counter terms Z = l+gn(n— l)/e, 
Z 3 = 1 + g R (7 - 4n)/e and Z x = 1 + 3g(n - l)/4e. They 
do not depend on / at any order since they correspond 
to "bulk" divergences for the random walk in interaction 
with the potential at the origin ("boundary" term). For 
the counter term Z/ we need Si , whose Feynman diagrams 
were computed in eqns. (8)-(ll): 



Si = nLfe~ 



1 + 



g(n - 1) (n + 2 
2~T~ 



2n{q') 2 L 



2/ 
e 



(14) 



We absorb the poles by the counter terms of the force- 
free theory and Z/ = 1 + 2//e. This expression depends 
neither on the number of replicas n, nor on g. Thus at first 
order there is no coupling between force and disorder. A 
straightforward calculation gives the RG /3-functions 



0f(9R,fR) = 
9 (9rJr) = 



dlogp 
dgn 



+ /I) 



(15) 



dlogp 



= -eg R + (5 - 2n)g% (16) 



In the physical case of random RNA, n = 0, they yield 
the RG flow depicted in fig. 3 with four fixed points. The 
attractive Gaussian fixed point O = (0,0) describes the 
molten phase of [4]. D = (e/5,0) is the fixed point of the 
glass transition found in [9,10]. We identify F = (0, e/2) 
with the denaturation transition of a homopolymer (wet- 
ting in 1 + 1 dimensions). More interesting, a new bi- 
critical UV unstable fixed point B = (e/2,e/5) emerges 
from our one- loop analysis. It leads to four phases sepa- 
rated by the critical lines / = e/2 and g = e/5. In par- 
ticular, a new phase at both high force and strong disor- 
der emerges. Physically, it corresponds to isolated frozen 
branched structures separated by free parts of the strand. 
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Fig. 3: Schematic view of renormalization group flow in the 
(g, f) plane. The critical lines distinguish the molten phase 
(I), a glass phase (II), (III) the denatured homopolymer phase. 
A new phase (IV) corresponds to RNA strands with isolated 
frozen structures separated by free parts. 



The extension £ as a function of the force / is given by 
the anomalous dimension of this force. The scaling ansatz 
reads 

C{f) = L<*F{[f-f c ]V) , (17) 

where J-(z) is a scaling function. At the bi-critical fixed 
point, we find 



7 = - dim[/ - f c ] = 



df R 



fc 



(18) 



Since C = (^f^ dx A(x)^, the exponent a = 7. Demand- 
ing that for large systems the extension C become exten- 
sive, i.e. ~ L, yields T{z) <~ z~~ x for z > 1, so that for 
large RNA molecules £(/) - L(f - / c )^ _1 . Setting e = 1 
yields the exponent 7=1/2 which is the same result as for 
the homopolymer denaturation transition, first discussed 
in [11]. This is consistent with the fact that the 1-loop 
beta function for / does not depend on g, so that 7 is not 
changed by the quenched disorder. However, this result 
does not hold at higher orders [17]. 

Conclusion. — To summarize, we have developed a 
field-theoretic description of random RNA under applica- 
tion of an external force, by extending the force-free field 
theory. It permits to study the second-order denatura- 
tion transition. At 1-loop order, the RG flow functions for 
force / and disorder strength g decouple so that the denat- 
uration is not influenced by disorder. We have computed 
the critical exponent 7=1/2 for the force-extension char- 
acteristic at the transition which agrees with previously 
considered homopolymer models. 

We have extended our calculations to second order 
in perturbation theory [17]. The renormalization group 
flow of / then depends on the disorder strength g. The 



procedure yields a scaling exponent 7 = 0.6 for the 
force-extension characteristic, resulting in C(f) ~ L(f — 

We invoke the locking hypothesis [9] in order to con- 
jecture the critical exponents in the glass phase. The ar- 
gument relies on the exact inequality dim[^>] > dim[&\ 
for the dimensions of contact and overlap operators. For 
e = 1 consistency of the theory requires that the inequal- 
ity is satisfied. Physically this means that different repli- 
cas follow the same minimal-energy path already at the 
transition. This is quite unusual. Normally one expects 
several minimal-energy paths to exist at the transition, 
and the strong-coupling phase to be characterized by a 
collapse of these distinct paths into a single one, resulting 
into different physical exponents. In a situation where at 
the transition only a single minimal-energy path exists, we 
cannot have a collapse of several paths, and the accompa- 
nying change of critical exponents. In the absence of a 
different mechanism, the exponents will not change for a 
small increase in disorder, and by renormalization argu- 
ments in the whole strong-coupling phase. We conjecture 
that this is what happens for the RNA freezing transition. 
It is then tempting to suppose this hypothesis to hold even 
in the presence of an external force, i.e. on the critical line 
beyond the bi-critical fixed point. This assumption leads 
us to extend our prediction a = 7 = 0.6 to the glass phase. 
Indeed, this prediction proves to be in reasonable agree- 
ment with numerical simulations of Miiller et al. [11, 18]. 
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